A Critical Analysis of Goldberger-Wise Stabilization in Randall- Sundrum Scenario 
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The Goldberger-Wise mechanism of stabilizing modulus in the Randall-Sundrum braneworld,by 

■ introducing a bulk scalar field with quartic interaction terms localized at the 3-branes has been 
extremely popular as a stabilizing mechanism when the back-reaction of the scalar field on the 
geometry is negligibly small. In this note we re-examine the mechanism by an exact analysis 
without resorting to the approximations adopted by Goldberger and Wise. An exact calculation of 

^— i ' the stabilization condition indicates the existence of closely spaced minimum and a maximum for 

^> , the potential and also brings out some new features involved in the context of stabilization of such 

■ braneworld models. 

The gauge hierarchy problem, which refers to the vast disparity between the weak scale and the Planck scale, leads to 
, the well-known naturalness problem in connection with the Higgs mass in the standard model of elementary particles. 
\Q ' One can avoid such problem by invoking supersymmetry at the expense of including a large number of particles 
£SJ , (namely superpartners) in the theory, all of which so far unseen. 

t-H ■ Recently, braneworld models with large extra spatial dimensions have been suggested [l[ which show that a large 
extra dimensional volume may bring down the fundamental mass scale of the theory to the TeV scale and thus resolves 
the naturalness problem. In an alternative approach Randall and Sundrum (RS) Q proposed a higher dimensional 
scenario consisting of a small extra spatial dimension orbifolded as S1/Z2, with two opposite tension branes located 

■ at the two orbifold fixed points. With a negative five-dimensional bulk cosmological constant, they obtained the 
"^T* ' following space-time metric as a solution of the five dimensional Einstein's equation, 

■ ds 2 = e- 2krW r 1 ^dx ll dx v - r 2 d<t> 2 (1) 

Such a warped geometry immediately resolves the naturalness problem by warping the scalar masses from Planck 
scale in one brane (called Planck brane or hidden brane) to TeV scale in the other brane( called the visible 
brane or standard model brane) for very tiny value (~ near Planck length) of the modulus r which measures 
the separation of the two branes. In this scenario, r is associated with the vacuum expectation value (VEV) of 
a massless four-dimensional scalar field (modulus field) which has zero potential, so that r is not determined by 
the dynamics of the model. Goldberger and Wise (GW) [|[ proposed to generate such a potential classically by 
introducing a bulk massive scalar field with quartic interaction terms localized at the two 3-branes and finally 
obtained a value kr ~ 12 by minimizing the potential, without any fine-tuning of parameters. In this analysis 
the back-reaction of the scalar field on the background geometry was neglected. A similar analysis by taking 
full back-reaction into consideration have been studied where the nature of the warped geometry as well as 
the possible stability of the modulus have been readdressedQ- In the present note however we stick to the 
original GW approach where the back-reaction can indeed be ignored and re-examine the modulus stabilization 
mechanism more critically. We extend the analysis beyond the approximations made in Q and bring out the 
dependence of the stability condition on the values of the various parameters. We show the existence of a maximum 
along with a minimum of the modulus potential for appropriate values of the parameters leading to possible instability. 
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We begin with the Goldberger-Wise analysis, which specifies an action of the form: 

S = SG rav ity + S v i s + Shid ~\~ S sca l ar , (2a) 



where, 

S G ravity 



J d A x r d(f>VG [2M 3 R + A] (2b) 

S vis = J d A x^-g s [L s - V s ] (2c) 

Shid, = J d 4 Xy/-g p [L p - V p ] (2d) 

Scalar = \ J <?x J* VG(G ab d A <Z>d B $ - to 2 $ 2 ) - J ^g~ p \ p {^> 2 - v p ) 2 - J d^x^g^X^ 2 - v s ) 2 . (2e) 

Here A is the five dimensional cosmological constant, V^, V p are the visible and hidden brane tensions. 

<I> develops a (^-dependent vacuum expectation value, which has to be determined from the classical equation of 

motion, 

<9 (e- 4CT <9 $) = mVe- 4(T $ + 4- 4<T A t ,r$($ 2 - v 2 s )5{<j) - tt) + 4 e - 4CT A p r$($ 2 - v 2 p )6(<j)) (3) 

where a = kr\(j>\. In the bulk (where the delta functions are not relevant), the equation has a general solution: 

$(</>)= e 2G \Ae va + Be- ua ] (4) 



where v = ^/4 + m 2 /k 2 . The solution is now plugged back into the original scalar field action and integrated over cj> 
to obtain an effective 4-dimensional potential for r, of the form: 

V»(r) - k(v + 2)A 2 {e 2vk ™ - 1) + k{y - 2)B 2 {1 - e ~ 2vk ™) + \ s e' ikr7r [^{nf - v 2 } 2 + A p [$(0) 2 - v 2 p ]. (5) 

The coefficients should be determined by matching the delta- function terms at the boundaries, as they appear in 
Eq.Q. The resultant condition on A and B is obtained as follows: 

fc[(2 + v)A + (2 - v)B] - 2A P $(0)[$(0) 2 - v 2 p ] = 0. (6a) 
ke 2kr7T [{2 + v)e vkr " K A + (2 - v)e~ vkrlx B\ + 2A s $(7r)[$(7r) 2 - v 2 } = 0. (6b) 

At this stage The authors of have argued that in the large A limit, one can choose $(0) = v p and <&(7r) = v s 
as the minimum energy configuration. However clearly with this choice both the Eq.([6]) can only be satisfied for a 
non-vanishing value of r only if A is infinitely large. 

Before re-examining the analysis of [3J] for arbitrary value of A ,we first calculate the first and second derivative of 
the potential to find the exact extremization condition for the modulus. It may be noted that such conditions were 
obtained in Q using some specific approximations in the large A limit. Our result here is exact and perfectly general 
for all values of A. A long but straightforward calculation yields, the first derivative for the potential (assuming 
explicit r dependence of A and B through the Eq.© as, 

V£(r) = -4kTr\ s e- 4kr7r {<5>{ir) 2 - v 2 ) 2 - Ak 2 n[(v + 2)e 2ukr * A 2 + (2 - v) e - 2vkrn B 2 + (4 - v 2 )AB\. (7) 

where prime denotes differentiation with respect to r. 

Using the extremization condition for the potential (V$(r) — 0) from the above equation we obtain a simple but 
exact form for the second derivative of the potential as 

V#(r) = 4fc 2 7n4(2 + v)AB' + [v - 2)BA'} (8) 

The sign of R.H.S. of the above equation determines exactly whether the stationary value of the modulus r is a stable 
value or not. 

Assume that for arbitrary value of A ( not infinity) the boundary value of the scalar field at the two orbifold fixed 
points are <&(</> = 0) = Q p (r) and $(^> = n) — Q s {r). Now we can easily express the undetermined constants A and 
B in terms of the quantities, Q p and Q s as, 

A = Qs(r)e- 2 ° - Q p (r)e-™ 
2sinh(i/<7) 

b = Q^-QM'- 2 * (10) 

2 sinh(iver) 
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Substituting the expressions for A and B in Eq[S] we obtain, 



sinh(j/ a) 



2v 



2v ' Q s 



sinh(i/ a) 



Qp r 2 + V c (v-2)o | ^ ~ 2 „-(V+2)cr 



2^ 



_ 2Ap Qp (n 2 _ 2\ 



(11) 
(12) 



Further using Eg 11 21 and the expressions for A and B into equ[7] under stability condition we arrive at ( for X s ^ 0), 

fc 



where 



x = 



1 + 



KQ 2 S 



C z 



C 



Q s 2v ' 2v 
Now, it is easy to manipulate the expression given below from the Eq|13l that is 



2 + v (u-2)a , V- 2 -(u+2)a I n 
2v 2v ' 



kr 



1 



■In 



1 



tx{v - 2) 

for the stationarity condition i.e V$(r) = 0. Here 

C = 



Q P (r) 



(13) 



(14) 



(15) 



/ ^{(2 + ^)e 2 (^- 2 ) CT - e - 4ff (4 - i/ 2 ) + (2 - v)e~ 2 ^+ 2 ) a } 
y 1 {(v + 2)e( l '- 2 ) ff + (v- 2)e-(^+ 2 ) CT } 2 



(16) 



" The expression for kr in equ ll5|) . is an exact expression for 



It may be noted that in the large kr limit, C y . 2 

the stationary value of the modulus r and is valid for any value of the brane coupling constants A. 
Using the above results, the expression for V" becomes 

Akirve~ 2a 

F * W = ~ sinh(zAx) ^ piQ2p " + Xs(Q " ~ v2 s^Q's + ^KUQl ~ <){Q 2 p - v 2 p )Q p Q s ] (17) 

where 'prime'{'} denotes derivative with respect to r. 

At this point, we want to emphasize that so far no approximation has been made and all the results are exact. 

Stability Analysis 

We now re-examine the stability of the modulus r. 
Case I: A — > oo 

From the expression for the potential EqfSl the minimum energy configuration leads to the 



Qs — > v. 



s i ^°zp r u p 



(18) 



which are constants. This completely agrees withe result of [3j. However some interesting new features also show up 
in this analysis which we discuss now. 

In this limit the expressions for the stationary points become(in the large kr limit) 



kr 



1 



ir{v-2) 



In 



2v v p (r) 1 
2 + 1/ v„(r) 1 i l v-2 



v+2 



In the infinity limit of the coupling constant 
So, the Eg IPTl becomes 

^'(r)=-3 



J ) 



Q' p = 



sinh {yd) 



(19) 



(20) 



(21) 
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FIG. 1: F(kr,v,n,t,q) with kr taking the different set of 
values for the parameters t and q but keeping fixed values of 
v = 2.02 and n = 1.5 




X 



FIG. 2: Plot of F(kr, v, n, t, q) with kr taking the different set 
of values for the parameters n and v but keeping fixed values 
of t = 0.001 and q = 0.001 



Now, equ llll tells that in a wide range of parameter values, the expression \ p (Qp — v p ) is negative where as from 
the Eq[T21 A p (Qp — Vp) is either positive or negative according to the value of kr. However, for 



1 



hi 



2v v p (r) 



1 



v v s (r) 



(22) 



V"(r) > 0, that means kr+ is a stable point for the minimum of the potential. So, for the other value of fer_, the 
potential will be maximum. Clearly no extreme fine tuning of the parameters is required to get the right magnitude 
for kr . Here, one important point is that the expression for the value of kr is different from that of [3J] . This is why 
if we calculate the value of kr corresponding to the values of parameters m/k = 0.2 and v p /v s = 1.5 given in Q we 
get kr — 10.846 which is lesser than what authors have predicted. But, for instance, taking v p /v s = 2.3 and v = 2.02 
yields kr — 12.2504 for the minimum of the potential. The value of kr — 14.4993 corresponds to maximum of the 
potential. Clearly these two values are very close to each other. 



Case II: A is finite but very large, 

From the Eq[6j it is clear that if the brane coupling constant is large but finite, the value of the scalar field Q p on 
the Planck brane is lower than v p and approaches v p as the value of the brane coupling constant tends to infinity. 
But for the standard model brane it is just the opposite where the scalar field value on the brane Q s , approaches v s 
from above as the brane coupling constant approaches infinity from a finite value. This prompts us to calculate the 
corrections to the boundary scalar field values in the leading 1/A order correction. These become 



Q P (r) 



Q s (r) =v s + 



kv. 



-2a 



X p v p 4sinh (ya) 



kVn 



■2a 



\ s Vs 4sinh {va) 



V s 


' 2 + \(2-»)« 


v-2 


_V p 


[ 2v 


f 2, 


?E - « 




v-2 


Vs 




2v 6 



-(f+2)<7 



(23) 



Now, we are interested in studying the modified expression for Eq 1 151 under large but finite values of brane coupling 
constants upto first order in 1/A correction. A straightforward calculation using Eqsl23lleads to the value, (in the 
large kr limit) 



kr 



tt{v-2) 



In 




-is)kr7T 



V t 



-qn n 



2v 



(24) 
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where, v,n = V p /v s , t — k/(X p Vp) and q — k/(X s v^) are four parameters of the model. Here also as in the previous case 
we obtain one minimum as well as one maximum for the potential. The positive sign corresponds to the value of kr 
which minimizes the radion potential. For this positive sign, the Eq. 1241 reduces to F(£x, v, n, t, q) = 0, where F is an 
appropriate function of the parameters. 

We plot the function F with respect to kr for different values of the parameters n, is, t, q as shown in the figures fT|2l 
and find out the roots of F(fcr, n, t, q) in the region where the value of kr solves the hierarchy problem i.e. kr ~ 12. 
The existence of nontrivial roots in this region is clearly shown in the above figures. 
Our analysis therefore brings out the following new features in the context of modulus stabilization: 
While Godberger-Wise calculation calculation was strictly correct only for infinite value of A our analysis here shows 
that this mechanism stabilizes the modulus r even for finite but large value of A. The stable value of the modulus in 
this case once again solves the hierarchy problem without any unnatural fine tuning of the parameters although the 
stable value for the modulus differs marginally from that estimated by GW. We have determined the modified value 
of kr , inflicted from the fmitcncss of A ,to the leading order correction around the infinite value i.e. in terms of inverse 
of A. As stated earlier our exact calculation indicates that there exists simultaneously a very closely spaced (~ Planck 
length(Zp)) maximum along with the minimum. This may have interesting consequences in a quantum mechanical 
version of such a model which however is beyond the scope of interest of this note. 
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